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Inductive Post Arrays in Rectangular Waveguide 

By T. A. ABELE 
(Manuscript received July 15, 1977) 

Previous attempts, based on mode -matching techniques, to obtain 
precise data for the equivalent circuit of inductive post arrays in rec- 
tangular waveguide have consistently failed due to convergence 
problems. A different formulation is presented for symmetrical post 
arrays, which is shown to be free from this defect. 



I. INTRODUCTION 

Waveguide band -pass structures employing cascades of inductive 
posts have been built for many years. They usually contain a symmetrical 
arrangement of posts in each cross-section, mostly one, two or three 
posts. The latter arrangement, for instance, is a favorite for X/4-coupled 
filters, since it strongly reduces higher order mode interaction. This al- 
lows A/4 spacings to be used instead of the 3X/4 spacings required for 
single post filters, thus leading to shorter filters. 

In the past all of these structures had to be designed on the basis of 
measured data for the equivalent circuit of the cross-sectional post ar- 
rangement, because the available theoretical calculations 1,2,8,9 are not 
sufficiently accurate. The obvious problem with measured data is, of 
course, that two errors are introduced, whose magnitudes are only poorly 
known: dimensional tolerances of the sample to be measured and errors 
in the measurement itself. 

Previous attempts to obtain theoretical data based on mode-matching 
techniques have consistently failed due to the convergence problem 
typically associated with taking a finite number of unknowns out of two 
sets of infinitely many unknowns. This paper presents a formulation 
which leads to only one set of infinitely many unknowns in the case of 
single or double posts. It may thus be expected that, when a finite 
number of these is taken, no convergence problem will be encountered. 
One may also speculate that this will continue to be true for arrays in- 
volving three or more posts, although in these cases more than one set 
of infinitely many unknowns is encountered again. 
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Fig. 1 — Post array. 

II. CONFIGURATION 

We wish to determine the equivalent circuit of the array in Fig. 1 in 
the plane z = 0. 

The posts are circular. They are numbered consecutively from \i = — M 
ton = M with n = designating the center post. The array is symmetrical 
with respect to the plane 2 = and the plane x = 0. The center post may 
or may not be present. Each post n has a diameter d M and a coordinate 
x = p M of its axis. Only dominant (TEi ) mode propagation is assumed. 
The surfaces shall be perfectly conducting. 

The electric field will be calculated as the superposition of two fields; 
the field which exists without the posts, the unperturbed field, and the 
field generated by the surface currents on the posts, the perturbation 
field. The surface currents, or rather the coefficients of their Fourier 
series, are treated as unknowns, which are subsequently determined in 
such a way that the tangential electric field vanishes on the surface of 
the posts. As usual, only two special cases of excitation are studied, even 
and odd, since this suffices to determine the equivalent circuit. 



III. UNPERTURBED FIELD 

We set 



with 



E 



TTX 



y even 



= (e>0* z + e "■"***) cos — 



TTX 



E y odd = (e jp ' z ~ e-JP* z ) cos — 



-W^-$ 



(la) 
(lb) 

(2) 
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Fig. 2— Post surface. 



where (dominant mode assumption) 

P—T >- (3) 

A a 

X is the wavelength in our medium. Obviously the fields in Eqs. (1) fulfill 
the boundary conditions everywhere except on the post surfaces. 

For later use we wish to develop these fields into Fourier series on the 
surface of a post v located at p„ and of diameter d v . 

From Fig. 2 we see that the post surface has the coordinates 



x = p„ - - d„ cos 



z = - d v sin 
2 



(4a) 
(4b) 



Introduction of these expressions into Eqs. (1) and use of the well-known 
expansion of e> 2 aine (Ref. 5, p. 22), results in 

£y even = L J n (\ /H,) c> n * 

X [ C08 (a P " ~" n *°) +(_1)ncos (~P" + n ^o)l (5a) 
E yoAd = t Jn(\pd v )ei"* 

n = — <x> \z / 

X l C ° S (a P " ~ n0 °) " (_1) " COS (a P " + " 0O ) 1 (5b) 



where 



eJ<t>o 



-J(* + 'l) 



(6) 
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IV. PERTURBATION FIELD 

To determine the field generated by the current distributions on the 
posts we observe first that all of these currents are independent of y and 
in the direction of y. Secondly, the effect of the broad waveguide walls 
can be replaced, making use of the common imaging technique, by as- 
suming that all posts are infinitely long in both y directions, and, again, 
have current distributions which are independent of y and in the di- 
rection of y. Thirdly, by employing the same imaging technique once 
more, we can replace the effect of the narrow walls by periodically re- 
peating the array of infinitely long posts in both x directions with post 
locations and current distributions, which are consecutive mirror images 
of each other. To determine the perturbation field, we can then simply 
sum up the fields generated by these infinitely many and infinitely long 
posts, without having to worry about the boundary conditions on the 
waveguide walls, since they are automatically fulfilled. 

From basic electromagnetic theory we get for the electric field gen- 
erated by a current filament stretching in y direction from -« to °°, lo- 
cated at x , 2 , and of strength I y , only the following component in y- 
direction 



joa/i r m e -./0l>/<*-*o) 2 +yo 2 +(«-zo) 2 | 
Ey = ~ "iT h J-|V(x-xo) 2 + Vo 2 +( 2 -2o) 2 l ^ 



OJjU 



= -^I y Ho™WW(x-x ) 2 + (2 -2o) 2 D (V) 

4 

The latter transformation may be found in Ref. 3, p. 27. n is the perme- 
ability of the medium. 

Making use of the symmetry of our structure and summing over all 
currents on all post surfaces we obtain 

•H- M -W - ( _ 1)W21 

4 m =-m Jo k=-°= 

(pW ( z --d M sin^) 2 + (x - p M + ka + ^d^cos^ |) d+ (8) 

Figure 3 explains the quantities 1^), d^, p M and the coordinates 

used. 

With this expression for the perturbation field we will do two things. 
First we will determine its value at a large distance to obtain expressions 
for the elements of the dominant-mode equivalent circuit. Secondly, we 
will evaluate it on the post surfaces in order to be able to come up with 
an expression for the boundary condition for the tangential electric field 
there. 
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Fig. 3— Post. 

V. PERTURBATION FIELD AT A LARGE DISTANCE 

If we write 

B——\z — cL sini£ 

a \ 2 * I 



C = a V ~ P " + 2 d " C0S ) 



A = 0a 

in Eq. (8), we can apply Eq. (34) from Appendix A to Eq. (8), which re- 
sults in 

M s*2ir 



u = -M JO 



E y = - 



i 



exp-;V^ 2 - 



(2/J - 1) 2 7T 2 



2 d„ sin \ff\ 

2 " 



r (2fe-i)x / i \-i 

X cos I x - p„ + - df, cos \p \ 



exp — 



k=i+i 



(2k - 1)V 



- R2 



z d^siiwf/ 



with 



1/ (26-1)^2 I 

fl |V a 2 -*| 

|~(2/z - 1)tt / 1 , \ 1] 

COS I a V ~ Pfl + 2 " C ° S V II* 



| 2 |>I dn , (2^111 </3< <2H^ 
la a 



(9) 
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The second of the two sums over y obviously represents the evanescent 
modes in the rectangular waveguide and, therefore, vanishes for large 
\z | . In accordance with our assumptions we have 

^</?<- (10) 

a a 

which means that / = 1 in Eq. (9). We therefore find from Eq. (9) for 

■KZ 

— »1 

a 

(3g(ln=-Mj0 

X cos I - (x — Pft + - dp cos \f/ j d\p 

= _ _wm_ g f 2 '/ M (iA) [e-#«*+M*-Pj.>/«e'i< w '' 8in ^^ 
2[i g a »=-m Jo 

+ e-i^gz-i^x-pjlaefydp sin (^-*o)] d\p (11) 

where (f> is again defined by Eq. (6). Using once more the expansion al- 
ready used in Eqs. (5) we obtain 

fi g a n=-Mm=-~ \Z ' 

X cos \- (x - Pfl ) + m0ol §* IM)e im *d* (12) 

The inversion of the order of integration and summation employed here 
presents no difficulty. 

Physical considerations tell us that the currents IJd) can be developed 
into a Fourier series. We write in the usual manner 

/,(*)- t c,^"* (13a) 

m= — " 

where, of course, 

c M . m = - L I hW"*** (13b) 

27T JO 

Eq. (13b), when combined with Eq. (12), results in 

_2^0M e _-, gZ M - ( _ 1)mcMmJfn /i^N 

X cob I - (x - p M ) - info (14) 
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Fig. 4 — Equivalent circuit. 

For x = (center of the guide) this results in 

(3 g a n=-Mm=-~ \2 v 



X cos (-p ll + m0o) (15) 



Combining this result with Eqs. (1) for x = we find for the total 
(unperturbed plus perturbation) field for irz/a » 1 

EySEF = eJpgZ * e~ jfig * ^^e~#** 

M ■ /l \ /7T \ 

x E £ (-i) m ^, mSd v 5 ^(-^) cos /- P/i + m0o ) (16 ) 

c ^.<J5 n ' s written here to distinguish between the values of c^ m for even 
and odd excitation. For the equivalent circuit of Fig. 4, which is valid for 
the plane z = 0, we obtain from Eq. (16) the reflection coefficient 

n even _ , i 2?raj P- 
Podd - ± 1 — — 

M /l \ /7T \ 

1 ^ ' - 1 ' ,: e lua sB r J m {- pd^J cos f - p M + mtf> J (17) 



V = —M m = — < 



with 



Peven = (18a) 

jX+ fs +1 

jX- 1 
Podd = - x+ (18b) 

These equations permit the calculation of X and R once the values of 
c /i,"iSdd n are known. 
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We note, that because of the structural symmetry with respect to 
x = 

C-,, m = (-D m C»,-m (19) 

for fi ■ 0, ±1, ±2 . . . Also, since we have 

/^«n(^) = ±/^en(27r-^) (20) 

it follows that 

Eqs. (19) and (21) permit reduction of c M , m for negative values of a* and/or 
m to those with positive values. 

VI. PERTURBATION FIELD ON POST SURFACES 

Referring once more to Fig. 2 we get for the perturbation field on the 
surface of a post v from Eqs. (8) and (4) 

/.in M /* 2tt ro . 

Ey=-^ L f hW £ <-U* 

y 4 „=-m Jo fc— • 

X H (2) (j8| (^ d, sin - - d M sin *) 

11 \ 2 l 1/2 \ 

+ (p v -p (1 + ka--d v cos<t> + -d> l cosyPj | j <ty (22) 

We wish to write for this a double Fourier series with and \J/ as inde- 
pendent variables. This can be done with the aid of the so-called "ad- 
dition" theorem (Ref. 5, p. 361) if we impose the condition that the posts 
do not penetrate or touch each other or the narrow walls of the wave- 
guide. We obtain 



MM M /»2ir 



P = 



x[£ £ £ (-D m+ *«/ m (£#,) j n (fa,) 

*a+p,-p M >0 

x ffJIL |0(p„ - P M + *o)l •#•*+"«« 

+ £ £ £ •^(J^)^^') 

fea-p„+p„>0 

X tf 2L |0(p M - p, + fca)} e ;(»*W) (-1)"+* 

+ Hm £ Jn(l0d v )H™(K{Sd v )e^*-»]d+ (23) 

«-.l„=-co \2 / \2 / J 

/i = »> only 
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Based on physical reasoning (summing contributions of current fila- 
ments in different order, integrating around each post before summing) 
we now exchange the order of summations and integration in Eq. (23) 
and carry out the integration. This leads to 

CO U IT CD /I \ l~ M oo /I \ 

E y J" E Jn(-(3d,)eJ»+\ L £ J m (-fid,) c„, m 

X £ (-l) k H%L n W\p„-p l , + ka\)[ S gn(p tl -p l , + ka)]"+'" 
k=— ™ 

+ 

With the abbreviation 



limH^^^Cn] (24) 



t m (-!)*«? (|A + *B|)Ngn (§ + *)]' 



A + kB *0 

= f m (A,B) = (-l)-/_ m (A,B) = (-l)™/ m (-AJB) (25) 
we can write this as 

y 2 



n t m Jn fad) e*+ 

[ m /l \ 

E E ^(-^)c,, m / m _ n {/3(p,-p 1 ,) I /?a} 

+ HmffJ» (^/?rf,)c^l (26) 

If we take this result for the perturbation field, add it to the incident 
field Eqs. (5) and impose the condition E y = on the surface of the post, 
we get (letting k — *■ 1) 

£ Jn (- fid) e>«4cos ^p„ - n<£ ) ± (-1)" cos /-p„ + n0 o ) 1 

X c MfTOar / m -„ |/8(p M - p„), 0a } + tf < 2) (| 0d,) c r ,„arl (27) 

Because of the uniqueness of Fourier series (Ref. 4, p. 186), this results 
in 

cos l-p ¥ - n0 o ) ± (-1)" cos l~p p + n<j>o) 

r M /i \ 

E E e7 m (-/SdJc M<maP / ra _ n { j 8(p M -p (r ), j 8a} 

+ «?(£*«,) c^ar] (28) 
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CL)/i7T 

2 



This equation holds for n = 0, 1, 2 It expresses the boundary con- 
dition on the surface of post v for even and odd excitation. If applied to 
all posts v ■ 0, ±1, ±2 . . . ±M, it expresses the boundary condition on 
all posts. However, because of the symmetry involved, only v = 0, 1, 
2. . M are needed. As before, Eqs. (19) and (21) permit reduction of c M>m 
for negative values of n and/or m to those with positive values. In sum- 
mary we can say that Eq. (28), if applied for n = 0,1,2,.. .and v = 0, 1, 
2. . M, will allow us to compute all of the unknown coefficients c„ >m . In 
turn, Eq. (17) will then allow us to compute the elements of the equiva- 
lent circuit, which means that our problem is solved. 

Appendix A and Appendix B provide expressions suitable for the 
computation of f m {A,B) in Eq. (28). These alternate expressions are 
essential, because the defining series Eq. (25) converges very slowly, as 
the magnitude of H% ] (z) decreases only with |z" 1/2 | for large z. The 
derivation of these expressions constitutes the most difficult and labo- 
rious part of this analysis. For convenience the results are repeated below 
in the form most appropriate for Eq. (28). From Eqs. (40) and Eq. 
(41), 

f m (pp, 0a) = - tan O cos (^ - ~) e^^o-w 
7r \ a I / 



O-tt/2) 



TT n =3,5. . . 



2\ /mrp rmr\ 

— cos ( — - J 

a \ a 2 / 



VfHE+lveH]' 



. 1 tr (m-n-l)\ 
• / 7r n = M!(^-2n-l)! 

m-2n 
A 



X 



. 7TP 

a sin — 
a 



h m -2n-i (cob^J (29a) 



f 2m -iWNa, 0a) = 



(29b) 



f 2 m(0Na, 0a) = - tan o (-l)^+^2m(« o -W2) 



+ ;i e 



a 



7r„=3,5..l //nA\2 J |~" X , L// nX \ 2 ill 2 " 

.1 ™ 2(-l) N+n (m + n- DK2 2 "- 1 - l)B 2n /X\ 2 " 



+J-E 

7T ra = 



(ro-n)!(2n)! 
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(29c) 



From Eq. (34e) and Eq. (34f) 

Trp 
cos — 

foWp, /3a) = eJ^-^/2) + j i l n ( cot 2 ?£) 

ir cos 0o 7r \ 2a / 



4 cos 



7T n = 3,5. . . 



rnrp 
a 



foWNa, /3a) = ii=^ e y(*o--/2) _ (_i)JV + yl / c + ^ 2a\ ( _ 1} 

7T COS 00 7T \ X / 



(29d) 



/V' 



+ /i E 



_11iV 



4(-l) 



(29e) 



'»--iver-i|[g + |v(S^i]" 

These expressions are valid if 

2a 
2a > X > — 

3 

p^O, ±a, ±2a . . . 

m = 1,2, 3, ... 

N = 0, ±1, ±2 . . . 

The polynomials /i m (x) are defined in Appendix B, Eq. (42). B n is the 
nth Bernoullian number and C is Euler's constant. 



VII. NUMERICAL RESULTS 

A series of calculations was made to investigate the question of con- 
vergence and to ascertain that the rather involved analysis is error-free. 
To this end the reactances X - 2/B and -X [Eqs. (18)] were calculated 
for the cases of single, double, and triple posts with X/a = 1.2, px/a = 0.25 
and 0do = /?di = 0.2 and 0.4, employing increasing numbers of variables 
and equations. Furthermore, the computed results were compared with 
measured data where such data were available. Lacking a full-fledged 
computer program the calculations were carried out with the aid of a 
programmable desk calculator, except for the matrix inversion, for which 
a general-purpose computer program was used. 

Table I summarizes the results of this work. The first observation that 
can be made is that, as expected, the convergence obtained for single and 
double posts is excellent. Three terms in the Fourier series for the post 
currents is all that is needed to obtain six place accuracy for the reac- 
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Table 



x B 



Number of posts 



a 



n. 



-X 



0.2 



0.4 



0.25 



0.25 



— 0.2 



0.4 



0.2 0.2 



0.4 0.4 





1 

2 

3 

4 
Meas. 



1 

2 

3 

4 
Meas. 


1 

2 
3 
4 

1 
2 
3 
4 



1 

2 

3 

4 
Meas. 



1 

2 

3 

4 
Meas. 



1.121835970 

1.121835438 

1.121835438 

1.12 

.6546985053 

.6546719818 



.6546719813 
.655 

1.121835969 
1.100680471 
1.100679708 
1.100679707 
1.100679707 
.6546985046 
.6071886432 
.6071426563 
.6071423836 
.6071423830 

.2599117670 

.2578995041 

.2578489230 

.2578488656 

.2578488652 

.255 

!02634303527 

.02329626388 

.02238859942 

.02238444183 

.02238437265 

.0265 



.009450748398 

.009450749381 

.010 

.03659710000 

.03659716655 

.037 

.009467748389 
.009467748828 
.009467799814 
.009467799806 

.03659710003 
.03685942441 
.03685949132 
.03685949138 

.01871192452 
.01872800993 
.01872801335 
.01872801340 
.020 

.07043107151 
.07063929772 
.07063993488 
.07063993658 
.074 



tances, which is more than enough for any technical application. The 
second observation is that, as was hoped, excellent convergence continues 
to exist for triple posts, even though in that case two sets of infinitely 
many unknowns are encountered instead of just one. Even for posts with 
susceptance values as high as B = 20 no more than four terms in the 
Fourier series are needed to obtain five-place accuracy. Presumably the 
analysis will converge even for four or more posts, but these arrange- 
ments are of little technical interest and thus probably not worth in- 
vestigating. Finally, when comparing the computed values with mea- 
sured data obtained with the aid of a very precise, computer-operated 
transmission measurement set, 10 sufficient agreement is found to as- 
certain that the analysis is free from any fundamental error. 



APPENDIX A 

We study the following series 



f(z, B, C,t) = E 



e -By/[(2n-l)*+t]*+z* CQS {[( 2n - 1)^ + t]C] 
V[{2n- l)lC+t\ 2 + Z* 



(30) 
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with t as a real variable, B > and real, C real, Rejz) > 0, Re 
|V[(2n - l)7r + t] 2 + z 2 \ > 0. This function is even in t and also periodic 
in t with the period 2ir. For reasons which will become apparent later 
we wish to develop it into a Fourier series in t 

~ao+ Y dkcoskt (31) 

2 k=i 

Without going into the fairly laborious detail the result is 



f(z,B,C,t) = - Y (-l) k K (z\VB 2 +(C + k) 2 \ cos kt (32) 

7Tfe=-co 

for the conditions stated for Eq. (30) plus either B > or B = and C 
y£ 0, ± 1, ± 2 Setting t = leads to 

- e -fiv/(2n-l)VM- z 2 cos ^ (2n _ 1)C j 

2 „?i V(2n - 1) 2 tt 2 + 2 2 

1 



= - E (-l)*K («|VB z +(C + fe) 2 |) < 33 > 

7T/ ? =- = 

provided Relz} > 0, Re(v / (2n - l) 2 x 2 + z 2 } > and either B > and real, 
C real or B = 0, C ^ 0, ±1, ±2 ... and real. The validity of Eq. (43) can 
be extended to include Re|z) = by analytic continuation. In doing this 
the points z = and z = ±j(2n — l)ir obviously have to be excluded, since 
at these points individual terms of the sums involved are not analytic. 
The result of the analytic continuation is for z —- jA 



Y. (~l) k Hg ] (A\VB 2 + (C + k) 2 \) 

fe=-« 

a, p-B-v/(2rc-l)27r2-A 2 

= 4j L = cos [(2n - DttC] (34a) 

n = iV(2n - l) 2 7r 2 - A 2 



with 



A > 0, real 
A ^ 7T, 3tt, 5tt . . . (34b) 



and either 



or 



argv/(2n - l) 2 *- 2 - A 2 = or - 



B > 0, real 
C real (34c) 

B = 
C * 0, ±1, ±2..., real (34d) 
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For the latter situation [Eq. (34d)] Eq. (34a) may be rewritten in the more 
rapidly converging form 

£ (-l) k H^(A\C + k\)=j-\n(cot 2 KA)+Aj £ 

/z = -a> 7T \ Z / n =l 

A 2 cos [(2n - 1)ttC] 

(2n - l)w V(2n - 1) 2 tt 2 - A 2 [(2n - 1)tt + V(2n - 1) 2 tt 2 - A 2 ] 

(34e) 
We also need a result corresponding to Eq. (34e) for C = N = 0, ±1, 
±2 . . . . We observe first that in the left hand sum the term k = —N has 
to be excluded for obvious reasons. Furthermore it is 

£ (-l)*Hl?\A\N + K\) = i-lV £ (-1)*HJ 2) (A|*|) 

= 2(-l)" £ (-l)W(^) (35) 
k=i 

An alternate expression for the last sum is known (Ref. 6, p. 333). We 
get 

.!/«., A\ 



k 

h * -N 



+4; n?i l_V(2/i - DV^A 2 ~ 2^rJ | 
= (-l)^(-l + 2;-(c + ln-)+4; £ 

I 1C \ 7T/ n =\ 



A 2 



(2n - l)wV(2n - 1) 2 tt 2 - A 2 [(2ai - 1)tt + V(2rc - 1) 2 tt 2 - A 2 \ 

(34f) 

with N = 0, ±1, ±2 . . . and Eqs. (34b) in force. Note that C in this last 
formula is Euler's constant. 

APPENDIX B 

We study for m = 1, 2 . . . the series 

fm(z,C,t) = 



z m cos ) \(2n - 1)tt + t\ C - m -J 



n=W[{2n - 1)tt + £] 2 + z 2 [(2n - 1)tt + t + V[(2n - l)ir + t] 2 + z 2 \ m 

(36) 
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with t as a real variable wi th the range ~w < t < -k, C real Re{z} > 0, 
Re{A/[(2n — l)7r + t] 2 + z 2 \ > 0. Analogous to the situation in Appendix 
A we wish to develop f m (z, C, t) + f m (z, C, —t) into a Fourier series 

— + 2_, a k cos kt 
2 k=\ 

over the range — -k < t < it, thereby continuing it periodically beyond that 
range. Again omitting the fairly laborious detail we get for C ^ 0, ±1, 
±2... 

f m (z,C,t)+f m (z,C,-t) 

= ^ (-l)™ +1 [4SmO'Cz)e' W2 - 8X m (|C| 2 )(sgn C)»] 

+ -H) m+1 £ (-l)*[2S m 0'C2+^2)e^-/ 2 

47T fe=l 

+ 2S m 0'Cz -jkz)eJ m "> 2 - 4K m (\C + k\z)\sgn (C + k)] m 

- 4K m (\C - k\z)\sgn (C - k)}™] cos kt = -±- (-1)™ 

Ztt 

X f {-l)k[2K m (\C + k\z)[sga(C + k)Y» 

- eJ mw/2 S m (JCz + jkz)] cos fet (37a) 
and for C = N = 0, ±1, ±2 . . . 

/ m ( 2) N,t)+/m(2,iV,-0 

= ^-(-D m £ (-l)*[2K m (|N + &|2)|sgn(iV + *)} m 

27T *?*-<» 

-eJ m * /2 S m (jNz + jkz)] cos fet + (-1) N cos — cos Nt (37b) 

m7r 2 

In these equations S m (z) denotes Schlaefli's polynomial (Ref. 5, p. 313). 
For t = this results in 



z m cos 



[~(2n-l)irC-ro-l 



„= i V(2n-l) 2 7r 2 + 2 2 [(2m - l)ir + V(2n - 1) 2 tt 2 + z 2 ] m 

= ^-(-l) m t (-l) k [2K m {\C + k\z)\sgn(C + k)Y» 

- e> m ' /2 S m {jCz + jkz)] (38a) 
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provided C^O, ±1, ±2 . . . , and in 

z m cos (2n - l)irN - m - 

? — 

»-i V(2n - l) 2 7r 2 + 2 2 [(2n - 1)tt + V(2n - l) 2 7r 2 + z 2 } m 

= y-(-D m £ (-l) fe [2^(|N + fe| 2 ){ S gn(N + ^)|- 
Air k=-<» 

- e jm-n S m (jNz + jkz)] + — (-1)" cos — (38b) 

2m?r 2 

for N m 0, ±1, ±2 .... By working on the two series with Schlaefli's 
polynomials the following alternative expressions are obtained 



z'" cos 



|(2n - l)xC-m-| 



n=\ V(2n - 1) 2 tt 2 + 2 2 [(2n - 1)tt + V(2n - 1) 2 tt 2 + z 2 ] m 

= ^-(-D m L (-l)*K m (|C + A|2)[8gn(C + A)]'» 
2ir * — a, 



+ J_ e ,W2 <g/ 2 (m-n-1)! /27TX m-2» 
4tt n =o n\(m — 2n — 1)! \jzl 

X I m g . . (39a) 

Ldx m_2n_1 Sin X Jx=Cr 



£ (-l)*X 2 x_ 1 (|N + fe|2)sgn(N + fe) = (39b) 

2 2X(_X)N+X 



2X 



„-! V(2n - 1) 2 tt 2 + z 2 [(2n - 1)tt + V(2n - 1) 2 tt 2 + z 2 ] 

-7" £ (-l) fc K 2X (|iV+fe|2) + -^-(-l)^^ 

27Tfc = -a> 4A7T 

. J-^x £ (Hfe-DK^-lfa ^ (39c) 
2x ' fcti (X-*)!(2*)! Vf/ 

where, again, C^O, ±1, ±2 . . , real, A/ = 0, ±1, ±2 . . . , m = 1, 2 . . . , 
X = 1, 2 ... , ReJ2} > 0, Re{V(2n - 1) 2 tt 2 + z 2 \ > 0. Following the same 
argumentation as in Appendix A the validity of Eqs. (39) can, by analytic 
continuation, be extended to include Re|2J = with the exception of 2 
= and 2 = ±j(2n — l)ir. The result is for 2 — *■ jA with A > and after 
some rearrangement 
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E (-l)*ffi 2) (|C + fc|A)[sgn(C + *)]'» 



k=-> 



-4/ £ 



A m cosT(2n -DttC-^1 



n =! V(2n - 1) 2 tt 2 - A 2 [(2n - 1)tt + V(2n - 1) 2 tt 2 - A 2 ] m 

7r n =o n!(n - 2n - 1)! \ A / 

X I ° 2 . . (40a) 

Ldx m-2n-1 sin x Jx=dr 

E (-l)*H^_ 1 (|N + fe|A)sgn(N + fe) = (40b) 

_E (-lW^(|N + /z|A) = 4y(-l)" + - 
XE 



1 2m 



„=i V(2n-l) 2 7r 2 -A 2 [(2/i - 1)tt + V(2n-l)V 2 -A z ] s 

+ i( -l)^ E (^ + ^-DK2^-l)B 2 „(-l)- /2^*» 
7V ' T n % (m-n)l(2n)\ \A/ 

valid for m = 1, 2 . . . , N = 0, ±1, ±2 . . . , C ^ 0, ±1, ±2 and real, A > 
0, real and A ^ ir, 3tt, 5tt . . . , and with arg \V(2n - l) 2 ir 2 -A 2 } = or 
7r/2. Eq. (40a) can be written in the following more convenient form 

E (-D fe //^ 2) (|C + / e |A)[sgn(C + *)]" 



k=-> 



A m cos (2n - 1)ttC - ro — 
= 4J „?, V(2n - D 2 7T 2 - A 2 [(2n - 1)tt + V(2n - DV^A^]" 1 

+ ; — E -J 777 ( l ■ r ) fc«-2n-i(cos Ctt) (41) 

47r „=o n\{m - 2n — 1)! \A sin Gx/ 

subject to the same restrictions as those enumerated for Eqs. (40). The 
polynomials h n (u) appearing in this equation are defined by 

h n (u)= (-l) n [sin" +1 o:-^--^— 1 (42) 

|_ dx n sin x Jcos x =u 

and can be shown to satisfy the following recursion formula 
which begins with 

h n (u) = nuh n -i(u) + (1 - u 2 ) dhn - liu) (43a) 

au 

h {u) = 1 (43b) 
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It appears impossible to give a closed form expression for these poly- 
nomials, but their coefficients can easily be calculated by the following 
scheme which is a consequence of the recursion formula 




6 61 



179' 



i.e., it is 

ho(u) = 1 

h\{u) = u 

h 2 (u) = l + u 2 

hs(u) = 5u + u 3 

h 4 (u) = 5 + 18u 2 + u 4 
It can be shown, incidentally, that the sum of all coefficients of h n (u) 
is equal to n! and that the coefficients of u° and u l are Euler's numbers. 
It furthermore appears, but has not been proven, that the coefficients 
are equal to those in Table 7.2.2 of Ref. 7, p. 260, the "number of per- 
mutations of the first N natural numbers with t u runs up." 
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